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The Apparent Size of a Closed Curve. 

By Arthur C. Lunn. 
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1. The apparent size of a closed curve in space, or the solid angle subtended 
by it as seen from any point, is defined as the solid angle of a cone with vertex 
at the point of sight and with the curve as directrix. This is measured, according 
to the definition of Gauss, by the area cut by the cone from a sphere of unit 
radius with center at the given point. It is equivalent to the apparent size of 
any surface bounded by the curve, provided a suitable convention be adhered to 
as to the algebraic sign of each element, according to whether one side or the 
other is visible from the point. Its value is then given by the familiar integral 

/cos (n, r) d S _ Pxcos(n, x) +y cos (n,y) + z cos (n t z)-, „ , >. 
— ~z -j ^ ^ i 1 ; 

extended over a surface bounded by the curve, with positive side at each point 
fixed by the direction of the normal n. Here the point of sight is taken as 
origin of coordinates and an element of solid angle is considered positive when 
the origin is on the positive side of the corresponding element of surface, 
r being the radius vector from the origin to a point of the surface. The explicit 
computation of o is required, for example, in the determination of the magnetic 
potential of a closed linear electric current.* 

The integral in (1) may be evaluated as the sum of three double integrals, 
extended over the areas bounded by the projections of the curve on the coordinate 
planes, or as one double integral extended over a plane area into which the sur- 
face is mapped by means of a representation of the coordinates of its points in 
terms of two independent parameters. But inasmuch as the surface here is 
purely auxiliary and is arbitrary except in that it is to be bounded by the given 
curve, it seems more natural to seek a method of computation such as to use the 
curve itself as realm of integration, for instance by integration with respect to a 
parameter t in terms of which the coordinates of the points of the curve may 
be expressed. 

* Gauss: " Allgemeine Theorie des Erdmagnetismus," Gesammelte Werhe, V, p. 170. 
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One such method is given by the ordinary formula 

o=/(l — cos$)dQ=f{\ — cos$)6'dt, (2) 

where 0, $ are the longitude and colatitude of the point where the unit sphere 
is intersected by the vector r from the origin to a point of the curve. Here the 
use of t as a variable of integration, in terms of which the coordinates may be 
considered as single-valued, relieves the ambiguity as to the interval of variation 
of 6 arising from the fact that the polar axis, being arbitrary, may or may not 
be linked with the curve. For comparison with other formulas to be considered 
it will be convenient to throw this into a form containing rectangular coordinates, 
which will here be done by means of a vector formula independent of any par- 
ticular coordinate system. 

Let a be a constant vector, of length a, in the direction of the polar axis 
implied in equation (2), and let Tj, ai be unit vectors in the directions of r and a 
respectively, so that T = rT 1) a = aa^ Then if t be thought of as time, 6' is the 
angular velocity of the plane of the vectors r, a about the fixed vector a as axis 5 
and cos $ = Ti . ai . * Now a vector normal to the plane in question is 
n = n Hi = T X a, so that the derivative ni of the unit normal vector is 
numerically equal to 6'. From this come: 

n' = rcn{ + n'iii = r'x a, Ji ! = (rxa) ! , nn' = (r' x a).(r x a), 
w 2 ni = n (r' x a) — Di {f x a) . (r x a), 
and from the last, by squaring and reducing : 

w 4 ni 2 = \{t x a) x (r' x a)f 8 = {[r'ra] a} 2 , (3) 

which shows that the vector of angular velocity of the plane may be taken as 

« = [ r ' ra l (_ a) (4) 

u (rxaf 1 aj ' w 

which is numerically equal to ni, but is parallel to the axis a while ni is 
perpendicular. This formula may be obtained also directly by a geometric 
construction, with r' treated as the velocity of a point moving along the curve, 
so that — [r'ra] /a is the moment of the velocity with respect to the axis. The 
sense of u, left undetermined by (3), is here so chosen that when r', r, a form 
a right-handed triple of vectors the direction of u shall be opposite to that of a. 
Thus if the forward circuit of the curve appear counter-clockwise as seen from 

* The vector notation is throughout that of Gibbs-Wilson : " Vector Analysis." 



188 



Lunn: The Apparent Size of a Closed Curve. 



what is taken in (1) as the positive side of the surface, the resulting value of the 
solid angle is 

J \ raj (r x a) 3 ' 
which in Cartesian rectangular coordinates is 



(5) 
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r = Vsc 8 + ?/ 3 + z 3 , 
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(6) 



- «3^) 2 + (« 3 « — «i 2) 2 + («i y — « 2 ») 8 
= a 3 r 8 — («i a; + a s y + a 3 z) 2 , 
where a 1; a 2 , a 3 are the scalar components of a. In this deduction the polar 
axis has any arbitrary direction with respect to the rectangular axes ; if it be 
taken to coincide with the z-axis, and t be the arc of the curve, this formula 
reduces to one used by Maxwell.* 

The integrand in this formula is independent of the particular set of axes 
used for the resolution of the vectors ; and the integral is independent of the 
special choice of the parameter t, but it contains the direction of the arbitrary 
vector a and thus depends apparently on the ratios of the scalar constants 
a i> a 2 , a 3 . The geometric source of the formula shows, however, that the 
integral around a closed curve must be independent of those constants, for any 
continuous variation of them which does not introduce singularities into the 
integrand. 

If the curve intersect the axial line, then at such a point r becomes parallel 
to a and the divisor n 3 vanishes. Since a is arbitrary, this could of course be 
avoided for any particular curve, but there is a feature of theoretical interest 
in such a case. 

If the intersection is on the positive half of the axis of a, so that $ = 0, 
the vanishing of other factors keeps the integrand finite; but if it is on the 
negative half, with ^ = n, the integrand becomes infinite. In the latter case 
a suitable value of a might be sought by a limiting process, for which purpose 
a short arc of the curve in the neighborhood of the intersection is replaced by a 
circular arc, for instance, with center on the axis. But the result then depends 
upon which way this arc lies with respect to the axis, there being in fact, for 



'Electricity and Magnetism," Vol. II, Ed. 3, p. 39. 
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arcs on opposite sides of the axis, two possible limiting values of the integral 
differing by 4n, which is the limiting value of the integral in (2) or (5) for a 
small circle with center on the negative part of the axis. Thus a variation of the 
vector a such as to make its negative continuation pass across the curve produces 
an abrupt change of 4n in the value of the integral. The obvious interpretation 
is that this 4n is the apparent size of the whole "sky," or the cyclic constant 
which occurs additively an arbitrary integral number of times in the general 
many-valued determination of a. Moreover, if two different expressions for the 
integrand in (6) be formed, corresponding to different choices of the vector a, 
their difference proves to be the exact derivative of a certain function of the 
coordinates which is expressible in terms of inverse trigonometric tangents 
and is many-valued with the same cyclic constant 4n. 

2. Another method of computation* depends on the relation « = In — a, 
generalized from a familiar theorem relative to polar spherical triangles, in which 
o is the length of that curve on the unit sphere which is polar to the intersection 
of sphere and cone. To translate this into an explicit general formula, let it be 
noticed that if the vector r goes to a point of the given curve then the unit vector 
of the vector p = r X r' goes to the corresponding point of the polar curve on 
the sphere. Then 

P p{ = r x r" -p' pi , / = (r x r') 2 , pp' = {vx r') . (r x r"), 

and from the first of these, by squaring and reducing: 

/pf = ?- 2 [r", r', r] 3 . (7) 

Here p{ is numerically the derivative of the arc of the polar curve, so that with 
a conventional choice of sign in taking the square root in (7) the solid angle is 

D = 27t— f r [*"'*' 'fi dt, (8) 

or in Cartesian coordinates : 



o = 2 n — f 



r 



x" y" z" 
x' y' z' 






dt. 



x y z 
p* = (y Z ' — Z y'f + (zx' — xz'f + (x y< — y x'f 

= r^q*— (xx' -f yy' + zz'f f 
2 3 = cc' 2 + y 2 + z' 2 . 

* Used lor example by Schwartz: Gesammelle Werke, I, p. 316. 
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To put this into a form using polar coordinates, for comparison with (2), 
let (x, y, z) be written (r sin $ cos 6, r sin $ sin 6, r cos $) ; then it reduces to 

n r($"d' — 0"A') sin <& — 20'd>' 2 cos$ — 0' 3 sin 2 $cos4> , , . 

The divisor p z vanishes at a point of the curve where the tangent line coin- 
cides with the radius vector. At such a point the integrand may remain finite 
if also the osculating plane pass through the origin ; but in general it becomes 
infinite and the value of a needs to be found by a limiting process, involving a 
system of curves obtained by warping the original so as to keepp different from 
zero. The typical case of such a critical point is that where the curve is such 
that its trace on the unit sphere presents a cusp. Then, according to whether 
the curve is deformed so as to change this cusp to a rounded arc or to some kind 
of loop, the limiting process gives a result indeterminate by an additive multiple 
of 47i, like that mentioned in §1. 

In contrast with the formula of § 1, the present one contains no arbitrary 
constants, but it involves the second derivatives of the coordinates and thus 
implies stronger conditions on the curve. For example, equation (5) or (6) 
could be used directly for a polygon or other curve having in general a definite 
tangent but with a suitably restricted set of points of discontinuity in the 
direction of that tangent; while in such a case equation (8) or (9) would be 
inapplicable unless for each such point there be added a corrective term 
measured by the corresponding change in direction of the line tangent to the 
trace on the unit sphere. 

3. The method now to be developed consists in a transformation of the 

surface integral of equation (1) into a line integral by means of Stokes's 

theorem : 

/curlv.da=/V.dr, (11) 

where dh is the vector equivalent of the scalar element of area dS. 

If R be put for the Newtonian vector —Ti/r 2 , then equation (1) may be 
written 

o=/R.tfa. (12) 

Thus if there is a vector function v such that curl v = R, then under the analytic 
conditions sufficient for Stokes's theorem the solid angle can be computed from 

u=fv.T'dt. (13) 
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Since R is a solenoidal vector function, there must exist for it such a vector 
potential V, whose determination in terms of rectangular coordinates involves 
the integration of the system of differential equations: 

dv 3 dv 2 x dv 1 dv 3 _ y dv % dv { z /, 

dy dz~~ r 3> dz dx~~ r H ' dx dy~ r 3 ' ^ ' 

which condition its scalar components (v lt v Z) v 3 ). 

By these equations alone V is not uniquely determined, but any two solutions 
differ by a potential vector V?7or one whose components are partial derivatives 
of a single scalar function U. This Uma,y be many-valued, but its discontinuities 
and those of V U must be allowed for by corresponding corrective terms in the 
application of Stokes's theorem. It may be predicted that the only indeter- 
minateness thus introduced will be such that the values of a computed by means 
of two different expressions for v will differ by a multiple of 4«. A general 
proof of this is not attempted here, but a verification is given below for some 
particular cases.* 

A particular solution of equations (14), found by a known tentative process, 
is the following : 
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(15) 



which may be written 

v = ir a xw; w = ( r 
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r + z 



)• 



(16) 



vr + x ' r + y' 

Then the general solution may be written v + VU, where the scalar function U 
is arbitrary. The special form of V just found gives, by substitution in (13) : 

« = i/[r', TuVlldt, (17) 

which in Cartesian coordinates is 
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r + x r -\-y r + z 



dt. 



(18) 



* For a theory of many-valued functions in a Riemann space see a paper by Dixon, Proceedings of the 
London Mathematical Society, Ser. 2, Vol. I, pp. 415-436. 
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With respect to this formula a qualification must be made, however, to cover 
cases where the vanishing of any of the quantities (r + x, r + y, r -f- 2) makes 
V infinite at any point involved in the integrals of equation (11). Since con- 
tinuity of V and of curl v is presupposed in Stokes's theorem, a corrective term 
is needed for every point where the negative part of a coordinate axis pierces 
the surface implied in (11), even when at all points of the curve itself the inte- 
grand in (18) remains finite. 

To find the magnitude of this correction it suffices to consider a special case, 
such as a circle with center on the negative 2-axis and lying in a normal plane, 

represented by 

x = p cos t, y = p sin t, z = — &, 

where p, b are positive constants, so that r is the constant Vp 2 -f b 2 . For the 
integral in (18) this gives 



1 
3 



_ r*' dt J p s bpsint bpcost 1 VI b\ 

r Jo "^ r — j r + p sin £ f + p cos t J * \ 3 r ) ' 



while the true value for a is 2n(l —bfr). Hence it may be inferred in general 
that for every point where the negative part of a coordinate axis, when taken 
in the direction from the origin to infinity, pierces the surface from its positive 
or negative side, a corresponding correction 47t/3, positive or negative respect- 
ively, must be added to the integral in (18). Thus the passage of the curve 
across a negative semi-axis during continuous deformation would be marked by 
a discontinuity of 4 71/ 3 in the value of the integral. 

For example, suppose o is the value of the integral for a curve linked 
with k of the negative semi-axes. By reversal of the sense of integration the 
same curve may be considered to define also the complementary solid angle 
with 3 — h corrective terms. These mutually complementary solid angles are 
then 6) + AnkjZ and — u + 4% (3 — &)/3, whose sum is 4 71, as it should be. 

If two curves lie on the same cone, with vertex at the point of sight, and 
are represented in terms of the same parameter t by a correspondence of those 
points which lie on the same linear element of the cone, then if the vector to 
one curve be r, the corresponding vector to the other curve is $r, where $ 
is a scalar function of t. On account of homogeneity of degree zero in the 
coordinates, the integrand in (17) is therefore the same for both curves, since 
the derivative 4/ has as coefficient a vanishing determinant. A similar remark 
applies to the alternative formulas (5) and (8). 
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Another solution of equations (14) is that obtained from (15) by reversing 
the signs of the coordinates, but not of r, and then changing the sign of the 
components of V to allow for the resulting interchange of right-handed and left- 
handed reckoning of rotations. The corresponding value of the solid angle in 
Cartesian form is 
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(19) 



which is to be subjected to a correction of 47t/3 for each piercing of the surface 
by a positive semi-axis. The value of V here used is obtainable from the former 
by addition of the vector whose components are partial derivatives of the 
scalar function 
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Another form is the mean of the two : 
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(20) 



if + s 3 z 3 + x 2 x 2 + y % 

which needs correction by 27t/3 for each piercing of the surface by either a 
positive or a negative semi-axis. This last form has the advantage that the 
irrationality r = \/ cc 3 + y % -\- z % does not occur in the determinant. 

The differential equations (14) have the same form for all systems of 
rectangular axes, but any particular solution for the vector v has components 
whose analytic form depends on the choice of axes. Thus if formula (18) be 
considered as written for one particular system of axes and then be transformed 
to another system, the resulting formula will contain the arbitrary constants 
which specify the relative orientation of the new system with respect to the old. 

If (x 1} y u %) be coordinates referred to a system of axes obtained from the 
system of (x, y, z) by rigid rotation, as represented by the equations 



X! = a n x + a 12 y + a 13 2, 
2/! = a 21 x -f a^y + a n z, 
z x =a 31 x + a 32 y + a 33 z, 



(21) 
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with the familiar conditions on the coefficients, then the integrand, formed 
according to the scheme of (18), but in terms of (x u y u z x ) and then expressed 
in terms of (x, y, z), is equivalent to 
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(22) 



where the second determinant factor is unity, and {x u y 1} z a ) are now to be 
understood merely as abbreviations. This may be obtained from the integrand 
explicitly given in (18) by addition of the exact derivative of the scalar function 
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Further alternative formulas for o could be obtained indefinitely by taking 
as integrand any arbitrarily weighted mean of expressions such as (22), each 
corresponding to a different transformation of axes, and demanding for the axes 
belonging to each term a similarly weighted correction for discontinuity. Of 
special interest among such would be those where the various systems of axes 
are related by means of the transformations of a certain rotation group, so that 
the resulting formula would be formally invariant under the transformations of 
that group, if each term of the integrand be given the same weight. 
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